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(54) A method for reconciling data at seismic and well-log scales in 3-D earth modelling 



(57) A geostatisical method for reconciling the dis- 
parity in scale between vertically detailed log measure- 
ments of a selected rock property in boreholes and 
vertically-averaged measurements of the same rock 



property as derived from seismic observations over a 
region of interest. 
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Description 

[0001 ] A statistical method for explicitly accounting for 
the difference in vertical scale between 3-D seismic 
data and stratigraphic measurements made using log- s 
ging tools. 

[0002] As is well known to geophysicists, a sound 
source, at or near the surface of the earth, is caused 
periodically to inject an acoustic wavef ield into the earth 
at each of a plurality of regularly-spaced survey sta- w 
tions. The wavefield radiates in all directions to insonify 
the subsurface earth formations whence it is reflected 
back to be received by seismic sensors (receivers) 
located at designated stations at or near the surface of 
the earth. The seismic sensors convert the mechanical is 
earth motions, due to the reflected wavefield, to electri- 
cal signals. The resulting electrical signals are transmit- 
ted over a signal-transmission link of any desired type, 
to instrumentation, usually digital, where the seismic 
data signals are archivally stored for later analytic 20 
processing. The travel-time lapse between the emission 
of a wavefield by a source and the reception of the 
resulting sequence of reflected wavef ields by a receiver, 
is a measure of the depths of the respective earth for- 
mations from which the wavefield was reflected. Certain 25 
attributes of the signals such as signal amplitude level 
and frequency may be related to the type of rock 
through which the wavefield propagated or from which 
the wavefield was reflected. 

[0003] The seismic survey stations of a 3-D survey are 30 
preferably distributed in a regular grid over an area of 
interest with inter-station spacings on the order of 25 
meters. The processed seismic data associated with a 
plurality of seismic traces from a plurality of receivers 
distributed over an area under survey may be formatted 35 
by well-known analytic methods to create a model of a 
volume of the earth. The model may be structural, 
showing the physical attitudes of the respective earth 
strata of an entire geologic sequence such as depth, dip 
and strike. Alternately the model may be designed to 40 
exhibit variations in the level of a selected seismic 
attribute, such as seismic impedance. Seismic imped- 
ance may be statistically related to variations in the tex- 
ture of some physical property of subsurface rocks such 
as porosity or lithology. The purpose of such earth mod- 45 
els is to select optimal locations for boreholes drilled to 
recover some desired natural resource such as fluid 
hydrocarbon products. 

[0004] A detailed picture of the stratigraphy of the ver- 
tical geologic column at a wellsite can be provided by so 
one or more logging tools of any one of many different 
known types, that may be lowered into a borehole on a 
wireline, on the drill stem in MWD operations or on 
coiled tubing. The vertical measurement resolution of 
wireline tools is very good; it may be on the order of cen- ss 
timeters or often, even millimeters. Horizontally, bore- 
holes are far apart, particularly in newly-developed 
prospects; the well spacing may be hundreds or more 
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meters. For that reason, the horizontal resolution of tog 
data measurements is very poor. 
[0005] In a 3-dimensional area! seismic survey, the 
seismic stations are closely-spaced horizontally, per- 
haps on the order of 25 meters or less, thus providing 
excellent lateral resolution. But the vertical, time-scale 
resolution of seismic data is a function of the frequency 
of the reflected seismic energy and the wavefield prop- 
agation velocity. The vertical resolution of typical low- 
frequency reflected seismic signals, after processing, 
may be on the order of tens of meters. 
[0006] It is useful to statistically relate the laterally 
sparsely-sampled but vertically detailed borehole meas- 
urements of selected rock characteristics with seismic 
observations which are densely spaced horizontally but 
poorly resolved vertically. In view of the limited vertical 
resolution of seismic measurements, seismic attributes 
are typically correlated with petrophysical well data 
averaged vertically across earth strata which may be 
several tens of meters thick. The seismic attributes are 
then used to guide the areal interpolation of the well- 
derived zone average data. The techniques of kriging 
and cokriging are commonly-used geostatistical meth- 
ods for performing this interpolation task. Kriging is a 
spatial prediction method for generating the best linear 
unbiased estimate of a rock property value by statistical 
interpolation of well data. Cokriging is a spatial predic- 
tion method using seismic attributes as secondary data 
to guide the interpolation of the average rock character- 
istic measured in wells. Both methods are analogous to 
Wiener filtering techniques applied in the time domain. 
[0007] The cokriging method is useful to derive an 
areal model of an average rock property for a given 
earth stratum by combining well and seismic data. How- 
ever, this and other similar techniques do not provide a 
vertically detailed 3-D model of the subsurface layer. 
There is a need for constructing such detailed 3-D mod- 
els from seismic and well log data. Three-dimensional 
modeling requires careful integration of log and seismic 
data at different scales. A review of prior integration 
efforts follows. 

[0008] S. B. Gorell, in a paper entitled Creating 3-D 
reservoir models using areal geostatistical techniques 
combined with vertical well data, 1995, SPE paper 
29670, presented at the Western Regional Meeting at 
Bakersfield, CA, and C. S. Burns et at. in a paper enti- 
tled Reservoir characterization by seismically con- 
strained stochastic simulation, 1993, SPE paper 25656, 
8th annual Middle East Oil show, Bahrain, propose an 
empirical scaling method. Here, each vertical column of 
cells in a 3-D porosity model is linearly re-scaled to 
reproduce a seismic derived average porosity map. The 
technique has the advantage of being straightforward to 
implement but the re-scaled 3-D model will not tie at 
deviated wells. Furthermore, the rescaling process may 
distort the data histogram. 

[0009] C. V. Deutsch et al, in Geostatistical reservoir 
modeling accounting for precision and scale of seismic 
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data, 1996, SPE paper 36497, Annual Technology Con- 
ference, Denver. CO, introduced an heuristic procedure 
based on simulated annealing. An objective function 
was constructed including a term that measures the 
degree of misfit between vertical average data and aver- 5 
age values computed from the 3-D model. Simulated 
annealing is used to perturb the 3-D model until the mis- 
fit is reduced to a value below a specified tolerance. The 
method has the advantage that is it possible to account 
for the precision of the average information. The method w 
is, however, very greedy of computer time and suffers 
from convergence problems when conflicting con- 
straints are present in the objective function. 
[0010] R. A. Behrens et al. in 1996, also at Denver, in 
SPE paper 36499 entitled Incorporating seismic 15 
attribute maps in 3-D reservoir models, teach a sequen- 
tial simulation technique based on conventional block 
kriging. A porosity estimate is obtained for each cell of a 
3-D model as a weighted linear combination of adjacent 
cell porosity data and the seismic derived porosity aver- 20 
age of the vertical column containing the cell. The 
method accounts explicitly for the difference of support 
volumes between well and seismic data but not for the 
fact that the seismic porosity averages may not be exact 
measurements. Compared to simulated annealing, this 25 
method has the advantage of being more robust and 
amenable to analytical analysis. But it involves solving 
relatively complex kriging systems constructed from 
average covariance functions. 

[0011] There is a need for a 3-D simulation method 30 
that will account for the difference in scale between 
seismic and logging data and which will be economical 
of computer processing time. This disclosure teaches a 
stochastic interpolation method which may be used to 
generate a vertically detailed 3-D earth model of a rock 35 
characteristic by combining finely-sampled log data with 
seismic attributes representing vertical average of the 
rock property of interest. 

[0012] A computer-aided geostatistical method for 
simulating, from seismic data, the spatial distribution of 40 
a designated rock property over a three-dimensional (3- 
D) array of cells. The array represents a discretized 
model of a subsurface earth layer wherein the vertical 
cell thickness is substantially less than the vertical reso- 
lution of the seismic data. A first set of log derived 45 
measurements of the designated rock property, made at 
control cells along well trajectories, is assembled in a 
programmed computer. A second set of measurements 
is assembled in said programmed computer corre- 
sponding to seismically derived, vertical average values so 
of the same rock property at each vertical column of the 
cells in the model. A cell to be simulated is randomly 
selected from the 3-D array. A search region is defined, 
centered around the cell to be simulated. The kriging 
estimate and variance of the rock property is calculated 55 
at the chosen cell from the rock property values at con- 
trol cells located within the search region. The most 
likely rock property value is determined at the cell based 
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on the measured vertical average for the column con- 
taining the cell and the rock property values at control 
cells in the same column. The variance of the Gaussian 
likelihood function is calculated at the cell under consid- 
eration. Thereafter, the mean of the Gaussian posterior 
probability distribution is determined for that cell as a 
weighted linear combination of the previously-derived 
kriging and maximum likelihood estimates. The vari- 
ance of the posterior distribution is now calculated from 
the kriging variance and the variance of the likelihood 
function. A simulated rock property is generated at the 
cell by sampling at random from the posterior probability 
distribution. The simulated rock property is entered into 
the array as an additional control cell. The entire proc- 
ess is repeated until all cells in the array have been sim- 
ulated. Thereupon, with the aid of a computer graphics 
program, the distribution of the simulated rock property 
is mapped over the 3-D array cells. 
[0013] The novel features which are believed to be 
characteristic of the invention, both as to organization 
and methods of operation, together with the objects and 
advantages thereof, will be better understood from the 



following detailed description and the drawings wherein 4££ . ^ 

the invention is illustrated by way of example for the pur- ^ ^ 
pose of illustration and description only and are not 
intended as a definition of the limits of the invention: v . . l ^ 

FIGURE 1 is a schematic vertical section through a , 

proposed 3-D earth model; 9 \ r . 

FIGURE 2 is a graphic illustration of equation (5); - $ 

FIGURE 3 is a flow diagram illustrating the best ^ 

method of operation. ; . ^ 

[0014] A discretized earth model is defined as a 3-D 



array of cells. Each cell is indexed by an integer, i. The 
vertical thickness of each cell is assumed to be much 
smaller than the vertical resolution of surface seismic 
data. Please refer now to Figure 1 which represents a 
vertical section through the discretized earth model, 
generally shown as 10. 

[0015] A rock property variable x s is defined in each 
cell, representing for example, porosity. After suitable 
transformation, the variables Xj are assumed to belong 
to a second-order, stationary multi-Gaussian random 
field with constant zero mean: 

E{x,}=0. (1) 

variance a 2 *, and spatial co-variance 

EiXjXjJ-Cu (2) 

where Cij is a function of the interdistance vector 
between cells / and 7. The cells in the earth model are 
assumed to have quasi point support, meaning that the 
co-variances in (2) are point-to-point values as taught 
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by to US patent 5,41 6,750, issued May 16, 1995 to P. M. 
Doyen et al., and which is incorporated herein by refer- 
ence. 

[0016] As shown in Figure 1 , let x 0 represent the rock 
property value to be simulated at current cell 0. Let z 5 
represent the seismic-derived vertical average rock 
property of the column of cells, 12, containing cell 0. 
The simulations must be constrained to reproduce the 
vertical average with a tolerance t; that is, 

10 

z = S j a j x j + t: (3) 
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in one cell is obtained by taking the product of the likeli- 
hood function, controlling the contribution of the seismic 
average information, and the conditional distribution of 
P( x ol x s). representing the influence of the previously 
simulated and original cell data. Compared to the likeli- 
hood function corresponding to collocated cokriging, it 
is worth emphasizing that the seismic average likeli- 
hood function introduced here depends on previously 
simulated rock properties in the column containing the 
currently simulated cell. 

It follows from our Gaussian assumption, that the first 
term in (5) is a Gaussian distribution given by 



where the sum is over alt cells in the column containing 
cell 0 and the Xj represent the simulated cell values. In 
(3), the coefficients aj represent user-specified averag- 
ing weights which may vary vertically from cell to cell as 
shown by graph 1 4 to the left side of Figure 1 . In (3), the 
error in each column, c, is modeled as a realization of a 
Gaussian white noise with zero mean and constant var- 
iance a 2 ,.. The variation of the seismically derived aver- 
age, z, from column to column is shown by graph 16 
above the earth model 10 in Figure 1. 
[0017] A simulated value for x 0 is obtained by sam- 
pling at random from the local posterior distribution 

p(x 0 |x s ,z) (4) 



15 p(x 0 |x s ) oc exp (6) 

I 2o SK J 



with mean x 0 equal to the simple kriging estimate of 
20 xq and variance ct 2 sk equal to the corresponding kriging 
estimation variance, the kriging being performed using 
Xg as conditioning data. 

[0019] In (3), the variable z was defined as the 
weighted sum of the variables in one column of the 
25 earth model. Ignoring the vertical correlation between 
the variables x in a column, it can be shown that the 
likelihood function l(x 0 ) is also a Gaussian kernel with 
mean given by 



where x s represents well measurements or previously 
simulated values at the set of cells as located within the 
search neighborhood 20 around randomly-selected cell 
0. In Figure 1, set s is equal to {1,2,3,4,5}. Sequential 
sampling of the conditional distribution (4) at all cells 
generates a simulated earth model constrained by the 
well data and the vertical averages. Ignoring terms 
which are not independent of x 0 and therefore not 
required in the sampling, (4) may be rewritten as 

p(x 0 |x s , z) cc p(x 0 |x s ) f(z|x c ,x 0 ) (5) 

where the Markov-type assumption has been made that 
f(z[x s , x 0 ) = f(z|x c ,x 0 ), 

with c the subset of s corresponding to conditioning 
cells in the vertical column containing cell 0. In Figure 1 , 
subset c={2,4}. This Markov assumption means that the 
conditional distribution of the column-average variable, 
z, depends only on cell values belonging to the column 
and not on the other cells. This is an approximation sim- 
ilar to that made in collocated cokriging. 
[001 8] When viewed as a function of x 0 for a given ver- 
tical average value z and fixed x c , t(x 0 ) = f(z|x c ,x 0 ) is 
defined as the likelihood function. The Bayesian updat- 
ing rule (5) merely states that the posterior distribution 
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35 and variance given by 

o^»(o^o| S a))/al (8) 

40 

• n ( 7 ). x o ML ma y be interpreted as the Maximum Likeli- 
hood estimate of x 0 based on the measured vertical col- 
umn average z and the weighted sum of all previously 
45 simulated values in the column xj, j c c. In (8), the sum- 
mation is over all cells not yet simulated in the column 
containing cell 0. 

[0020] Under the multi -Gaussian assumption, the 
local posterior distribution in (5) is also a normal distri- 

50 bution. This can be verified by observing that (5) is the 
product of two exponential kernels whose exponents 
are quadratic functions of xq. The product is therefore 
an exponential with an exponent quadratic in x 0 , defin- 
ing a Gaussian posterior distribution whose mean and 

55 variance can be shown to be 
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SK 2 ML 2 

X 0 a ML + X Q a SK 



(9) 
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SK 



and 
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y ML 



SK 2 
y ML 



(10) 
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A graphic illustration of the process for obtaining the 
posterior distribution as the product of the Gaussian 
likelihood function and the Gaussian, kriging-based dis- 
tribution is shown by curves 17, 18 and 19 in Figure 2. 
[0021] The best mode of using this invention, based 
on formulas (9) and (10) as presently contemplated is 
best shown by means of a flow diagram as illustrated in 
Figure 3: 

Step 1 . Assemble, in a programmed computer, log- 
derived data at control celts corresponding to well 
locations. 

Step 2. Assemble seismically-derived vertical aver- 
age data for all columns of the 3-D array of cells. 
Step 3. Select at random, a cell, 0, to be simulated. 
Step 4. Calculate x 0 SK and ct 2 sk by point kriging 
using the data Xg which correspond to all original 
and previously simulated data at control cells falling 
within a specified search region. 
Step 5. Calculate the most likely value x 0 ML and the 
variance of the Gaussian likelihood function a 2 ML , 
according to formulas (7) and (8). 
Step 6. Using equations (9) and (10) and the results 
of steps 4 and 5, calculate the mean m 0 and vari- 
ance a 0 2 of the Gaussian posterior distribution. 
Step 7. Obtain the simulated value x 0 by sampling 
at random from the Gaussian distribution N(m 0 , 
a 2 0 ) and consider the simulated cell value to be an 
additional control cell. 

Step 8. Iteratively return to step 3 until all cells in the 
3-D grid have been simulated. 
Step 9. When all cells have been simulated, map 
and display a tangible simulated global model with 
the aid of any well-known computer graphics pro- 
gram. 

[0022] Computer implementation of the simulation 
process with vertical average constraints as described 
above is particularly simple and efficient as it does not 
require any block kriging calculations. All that is required 
is a straightforward modification of the classical Gaus- 
sian simulation process to perform the Bayesian updat- 
ing of the kriging mean and variance according to 
formulas (9) and (10). The impact of the average data is 
easy to understand thanks to the decoupling of the influ- 
ences of point data and column average data. 
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[0023] Following are some important properties of the 
simulation procedure: 

[0024] The mean of the posterior distribution in (9) is 
a linear combination of point kriging estimate x 0 SK of the 
cell value and the maximum likelihood estimate x 0 ML 
based on the vertical column average and the previ- 
ously simulated cell values in the column. 
[0025] The simulations tie at the wells. This follows 
from the exactitude property of piping, and the fact that 
(9) and (10) reduce to m 0 = x 0 s and a 0 2 = 0 at cells 
corresponding to well data. 

[0026] The simulated earth models approximately 
reproduce the vertical column averages with a tolerance 
depending upon the value of the noise variance a, 2 . In 
particular, the vertical column average data are repro- 
duced exactly when the noise variance is zero. This 

2 2 

results from the fact that cr ML = a 0 =0 when the 
last cell in a column is simulated and the noise variance 
is zero. From (7) and (9), the last simulated value 
Xq is determined from the equation 
m 0 = x 0 ML = (z - Sa^xp/aQ , where the sum is over all 
cells in the column except the last to be simulated. In 
practice, exact reproduction of seismic-derived column 
averages is not desirable and may give rise, to abrupt 
changes in simulated values. A-laterally variable toler- 
ance level may even be specified to reflect changes in 
the accuracy of the seismic average information. 
[0027] The simulation procedure works for any linear 
average of the form given in equation (3). We may 
select equal weights a^ = 1/n in equation (3) for a sim- 
ple arithmetic average. On the other hand, y if seismic- 
derived vertical averages are more sensitive to rock 
property values near the top of a layer, we may choose 
weight values which decrease with increasing depth in 
the layer as shown by graph 14, Figure 1 . 
[0028] It should be observed that the simulation proc- 
ess reduces to point kriging simulation at vertical levels 
in the layer where the averaging weights are zero. Thus, 
when a 0 = 0, (9) and (10) reduce to m 0 =x 0 SK and 
a 0 = a SK respectively. 

[0029] In arriving at approximations (7) and (8), we 
have ignored vertical correlation in the calculation of the 
column average likelihood function. Exact expressions 
may be obtained by replacing the coefficients a in (7) 
and (8) by a+X, where the X are block kriging weights. 
When a i: is equal to zero and the vertical averaging 
weights aj are constant, the exact solutions are numeri- 
cally identical to conventional block kriging but are 
obtained in a simpler manner due to the decoupling 
between point data and vertical average data. 
[0030] In practice, experience shows that the approx- 
imations (7) and (8) may be applied without adversely 
affecting vertical continuity on the simulated earth mod- 
els. This is because the likelihood term is used to mod- 
ulate the kriging solution which already accounts for 
vertical continuity. 

[0031] This invention has been described with a cer- 
tain degree of specificity by way of example but not by 
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way of limitation. Those skilled in the art will devise obvi- 
ous variations to the examples given herein but which 
will fall within the scope and spirit of this invention which 
is limited only by the appended claims. 

5 

Claims 

1 . A computer-aided geostatistical method for simulat- 
ing, from seismic data, the spatial distribution of a 
designated continuous rock property over a three- w 
dimensional array of cells, the array representing a 
discretized model of a subsurface earth layer, the 
vertical thickness of said cells being substantially 
less than the vertical resolution of the seismic data, 
comprising: is 

a) assembling in a programmed computer 
including a data processor, a first set of log- 
derived measurements of the designated rock 
property made at control cells along well trajec- 20 
tories and a second set of measurements cor- 
responding to seismically-derived, vertical 
average values of the same rock property at 
each vertical column of cells; 

b) randomly selecting a cell to be simulated 25 
from the three-dimensional array and defining 

a search region around said cell; 

c) calculating the kriging estimate and the krig- 
ing variance of the rock property value at said 
cell from the rock property values at control 30 
cells located within the search region; 

d) determining the most likely rock property 
value at said cell based on the measured verti- 
cal average for the column containing said cell 
and the rock property values at control cells in 35 
the same column; 

e) calculating the variance of the Gaussian like- 
lihood function at said cell; 

f) calculating the mean of the Gaussian poste- 
rior probability distribution at said cell as a 40 
weighted linear combination of the kriging esti- 
mate and of the maximum likelihood estimate 
derived in steps c) and d) respectively; 

g) calculating the variance of the posterior dis- 
tribution from the kriging variance and the vari- 45 
ance of the likelihood function derived in steps 

c) and e) respectively; 

h) generating a simulated rock property value 
at said cell by sampling at random from said 
posterior probability distribution; so 

i) entering said simulated rock property value in 
said array to provide an additional control cell; 
j) repeating steps b) through i) until all cells in 
the array have been simulated; and 

k) with the aid of a computer graphics program, ss 
mapping the distribution of the simulated rock 
property values over the three-dimensional 
array of cells to provide a tangible global dis- 
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play of the rock property distribution. 

2. An article of manufacture comprising: 

a global display of a model of the spatial distri- 
bution of a simulated continuous rock property 
value over a region of interest as produced by 
the process of claim 1. 

3. A display system comprising : 

means for receiving input log-derived measure- 
ments of a rock property; and a set of measure- 
ments corresponding to seismically-derived, 
vertical average values of the same rock prop- 
erty; 

means for simulating the spatial distribution of 
the property using a method in accordance 
with claim 1 ; and 

means for displaying the property distribution. 



BNSOOCID: <EP 091 1652A2_I_> 




8 

BNSDOCID: <EP 091 1652A2J_> 



EP0 911 652 A2 




z 






o 






I- 






RIBU 




CM CO 

b 


h— 




+ 


CO 


_j 

CM 5 




a 


b 




z 




b 


< 






CO 






CO 


CO 


< 


b 


CD 






b 


o 








LU 












co 






o 












UL 






O 






LU 






O 






z 






< 


CO 




VAR 


b 

-J 

So 
X 




a 
z 


b 


< 


+ 


+ 


EAN 


b 


CMS 

b 








Ui 


coo 






X 




< 










) 


3 






(J 


o 


_J 


E 


< 






O 











a 




b 




O 






—i 






z 


LU 


LU 


a 


> 


— i 


h- 


o 


CO 


cr 
j— 


O 


z 


a 


o 




o 


O 


_l 


01 


< 


LL 


z 




o 


O 


h- 


Q 


Q 


Z 


Q 




< 










< 


— I 


X 


— I 
LU 


LU 


o 


VALU 


ATED 


Q 




LU 










CO 








§ 


CO 






h- 




Q 




Z 


a 


< 



Q O 




BNS0OCID: <EP 091 1652A2_I_> 




THISP4@IIUHK(usptb, 



(19) 



J 





(12) 



Europa^^Bs Patentamt 
European Patent Office 
Office europeen des brevets (1 1 ) EP0 911 652 A3 

EUROPEAN PATENT APPLICATION 



(88) 


Date of publication A3: 


(51) IntCI 7 : G01V 11/00 




30.08.2000 Bulletin 2000/35 




(43) 


Date of publication A2: 






OR riA 1QQQ Rulletin 1QQQ/17 
C.O*W** 1 DUIICllll 1999/1' 




(21) 


Application number: 98308059.9 




(22) 


Date of filing: 02.10.1998 




(84) 


Designated Contracting States: 


(72) Inventors: 




at rf r*v np nvc fq fh fr or rio if it i i i 1 1 

Ml DC Or! OT DC Li fx CO rl Fri OD urf IC 1 1 LI L.U 






MCNL PTSE 


Esher, Surrey KT10 8PU (GB) 




Designated Extension States: 


• Psaiia, David Edgar 




AL LT LV MK RO SI 


Ealing, London W1 3 9HS (GB) 


(30) 


Priority: 03.10.1997 US 943577 


(74) Representative: 






O'Connell, David Christopher 


(71) 


Applicant: 


Haseltine Lake & Co., 




WESTERN ATLAS INTERNATIONAL, INC. 


Imperial House, 




Houston, TX 77042-3115 (US) 


15-19 Kingsway 






London WC2B 6UD (GB) 



(54) A method for reconciling data at seismic and well-log scales in 3-D earth modelling 



CO 
< 
CM 

in 

CO 



(57) A geostatisicat method for reconciling the dis- 
parity in scale between vertically detailed log measure- 
ments of a selected rock property in boreholes and 
vertically-averaged measurements of the same rock 



property as derived from seismic observations over a 
region of interest. 



| SELECT CELL 0 TO SIMULATE AT RANDOM 



INPUT LOG-OERIVED DATA 
AT CONTROL CELLS 



CALCULATE KRIGING ESTIMATE 
AND KRIGING VARIABLE, 



INPUT SEISMIC DERIVED 
VERTICAL AVERAGE DATA 



CALCULATE MOST LIKELY VALUE AND 
VARIANCE OF LIKELIHOOD FUNCTION 

AND CTml 



CALCULATE MEAN AND VARIANCE OF POSTERIOR GAUSSIAN DISTRIBUTION 



mo * ■ 



-6 Fig. 3 



DRAW SIMULATED VALUE X AT RANDOM FROM POSITIVE N(m 0 .(To ) 
AND TREAT SIMULATED CELL AS ADDITIONAL CONTROL CELL r~ 7 




NO ^ ALL CELLS 
SIMULATED 



OUTPUT SIMULATEO MODEL AND DISPLAY 
USING COMPUTER GRAPHICS 



<3> 

O 

Q_ 
LU 



Primed by Xerox (UK) Business Services 
2.16.7 (HRS)/3 6 



BNSDOCID: <EP 091 1652A3_I_> 



EP0 911 652 A3 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 98 30 8059 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with indication, where appropriate, 
of relevant passages 



Relevant 
to claim 



CLASSIFICATION OF THE 
APPLICATION (lnt.Cl.6) 



US 5 638 269 A (FOURNIER FREDERIQUE ET 
AL) 10 June 1997 (1997-06-10) 

* column 2, line 9 - column 3, line 2 * 

* column 4 t line 62 - column 5, line 15 * 

US 5 416 750 A (DOYEN PHILIPPE M ET AL) 
16 May 1995 (1995-05-16) 

* column 2, line 55 - column 3, line 10 * 



US 5 321 613 A (VANDERGRAFT JAMES S 

AL) 14 June 1994 (1994-06-14) 

* column 3, line 56 - column 4, line 



ET 



62 * 



1-3 



1-3 



1-3 



G01V11/00 



TECHNICAL FIELDS 
SEARCHED (lnt.CI-6) 



G01V 



The present search report has been drawn up for all claims 



Place of saaicn 

THE HAGUE 



Data of completion of &* Match 

10 July 2000 



Examrwr 

Swart jes, H 



CATEGORY OF CITED DOCUMENTS 



X : particularly relevant if taken alone 

Y : particularly relevant it combined with another 

document of the same category 
A : technological background 
O : non -written disclosure 
P : intermediate document 



T : theory or pnncple underlying the invention 
E : earlier patent document, but published on. or 

after the ruing date 
D : document cited in the application 
L : document cited for other reasons 

& : member of the same patent family, corresponding 
document 



2 



BNSDOCID: <EP 091 1652 A3_r_> 



EP0 911 652 A3 



ANNEX TO THE EUROPEAN SEARCH REPORT 
ON EUROPEAN PATENT APPLICATION NO. 



EP 98 30 8059 



This annex lists the patent family members relating to the patent documents cited in the above-mentioned European search report. 
The members are as contained in the European Patent Office EDP file on 

The European Patent Office is in no way liable for these particulars which are merely given for the purpose of information. 

10-07-2000 



Patent document 
cited in search report 


Publication 
date 


Patent family 
member(s) 


Publication 
date 


US 5638269 


A 


10-06-1997 


FR 


2710418 


A 


31-03-1995 








CN 


1114511 


A,B 


03-01-1996 








DE 


69415142 


D 


21-01-1999 








DE 


69415142 


T 


06-05-1999 








EP 


0671017 


A 


13-09-1995 








WO 


9508781 


A 


30-03-1995 








JP 


8503783 


T 


23-04-1996 








NO 


952000 


A 


19-07-1995 


US 5416750 


A 


16-05-1995 


AU 


677757 


B 


01-05-1997 








AU 


1507195 


A 


05-10-1995 








CA 


2144347 A,C 


26-09-1995 








EP 


0674189 


A 


27-09-1995 








NO 


951029 


A 


26-09-1995 


US 5321613 


A 


14-06-1994 


NONE 









1 

o 

u. 

9 

uj For more details about this annex : see Official Journal of the European Patent Office. No. 12/82 



BNSDOCIO: <EP 091 1652 A3 J _> 



# 



THIS PAGE BLANK ww™' 



